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PREFACE 


This  report  documcnls  the  results  of  US  AFETAC  inlcrbranth  tasking  6075007»  Till  m  of  Houri>  Temperatures  for 
Long-Term  Missing  Data”.  It  describes  the  development  of  a  method  for  filling  in  missing  hourl>  temperature',  for 
kKdtions  that  report  only  daily  maximum  and  minimum  temperatures.  The  analyst  \vas  Mr  William  R  Sehaub,  Jr » 
USAFETAC/DNO. 

The  original  support  assistance  request  (USAFETAC  Project  60750,  "Request  for  Engineering  Weather  Daui 
Support")  was  directed  to  USAFETAC/ECE  from  HQ  USAF/LEEEU  for  the  calculation  of  c^Kjlmg  degree  hours 
tCDH).  Since  many  locations  that  require  CDH  calculations  have  incomplete  temperature  records,  ECE  asked  DNO 
to  provide  a  way  to  "fill  in"  missing  hourly  temperatures. 

This  study  is  an  extension  of  original  work  done  by  Maj  W.F.  Miller  (1990).  Miller  developed  a  method  for 
interpolating  hourly  temperatures  when  less  than  6  consecutive  hours  of  dau  arc  missing.  The  same  temperature 
data  was  used  in  both  studies. 
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1.  INTRODUCTION 


1.1  Purpose  of  the  Study.  USAFETAC’s  Engineering  Mct^^orology  Section  (USAFETAC/ECE)  asked 
USAFETAC/DNO  to  develop  a  method  for  estimating  missing  houriy  tenpcraturcs  for  locations  that  report  (ml> 
maximum  and  minimum  temperatures,  such  as  those  in  the  NO  A  A  Climatological  Data  (CD).  USAFETAC.'EC  E 
required  such  a  temperature  estimation  technique  before  they  could  calculate  cooling  degree  hours  (CDH)  for 
stations  at  which  only  maximum  and  minimum  temperatures  were  reported. 

1.2  Previous  Work.  This  study  is  an  ex.:insion  of  work  documented  by  USAFETAC/PR-90/()()6,  Shnn-lcrni 
Hourly  I cmpcraiurc  Interpolation,  by  Maj  Walter  F.  Miller,  December  I99().  In  his  study.  Miller  described  a 
method  for  interpolating  hourly  temperatures  when  less  than  6  consecutive  hours  of  data  were  missing  Input 
temperature  daUt,  the  .same  for  both  studies,  was  taken  from  stations  in  different  climates  to  test  the  temperature 
filling  methods  adequately.  Those  tests  arc  described  in  Section  3.  Test  results  and  a  description  ol  the  two 
computer  programs  developed  for  ECE  are  in  Section  4. 
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2.  TEMPERATURE  DATA  USED 


2.1  Database.  Hourly  .siirlaLC  icmpcniiurc  dal.i  ^d.s  obuiincd  liom  DB2  Uiblc.s.  The  viala  i.s  .stored  in  tenths  of 
dcftrecs  f  lor  worldwide  iniliUiry  <md  eivilian  rcporlii.i^;  stations.  Typical  period  of  record  (POR)  was  from  19.59  to 
19XK,  with  the  exception  of  .Scott  AFB,  which  only  had  hourly  temperature  records  from  1962  to  1983. 

2.2  Test  Stations.  As  U.SAFlirACVFCF.  recommended,  :hc  suitions  li.stcd  liclow  were  .selected  to  test  the 
lemperalure  lilliiit;  ukiIuhI.  as  sliowii,  they  arc  repre  .cniaiivc  oi  (.on’.incntal,  coastal,  dc.sert,  .subtropical,  and  tropical 
regions. 


STATION 

BLKSTN 

LAT 

LON 

ELEV  (Meters) 

REGION 

.Scot'  AFB,  11. 

724338 

38".33’N 

89".L3’W 

138 

Conlincnial 

Pease  /vFB,  Nil 

726055 

43"05’N 

70"49’W 

31 

Coniincnial/Coasial 

McCItord  AFB,  WA 

742060 

47"09’N 

22"29’W 

98 

Subtropical 

Luke  AFB,  AZ 

722785 

33"32’N 

112"23’W 

332 

Dcseri 

1  lomcsiead  AFB,  FL 

722026 

2.5"29’N 

90”23’W 

2 

Sublropical/Coasial 

Mildcnhali  AB,  UK 

03.5773 

.52"22’N 

0'’29’W 

10 

Continental 

Ramsiein  AB,  CiF. 

106140 

49"26’N 

7"  36’E 

238 

ContinenutI 

lionoiiiiii,  HI 

911820 

21”21’N 

1.57“  .56’ W 

5 

Tropical 

PERIODS  OF  RECORD 


.Scott  AFB 

Jan  1962  to 

Pease.  AFB 

Jan  1959  to 

McChord  AFB 

Jan  19.59  to 

Luke  AFB 

Jan  1958  to 

1  lome.siead  AFB 

Jan  19.59  to 

Mildcnhali  AB 

Jan  1959  to 

Ramsiein  AB 

Jan  19.59  to 

Honolulu 

Jan  19.59  to 

Aug  1985 
Dee  1988 

Dec  1970;  Jan  1973  to  May  1988 
Dec  1970;  Jan  197.3  to  May  1986 
Dec  1970;  Jan  1973  to  Dec  1988 
Dec  1970;  Jan  1973  to  May  1988 
Dec  1970;  Jan  1973  to  May  1988 
Jim  1988 
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3.  METHODOLOGY 


3.1  Previous  Models.  The  diurnal  variation  of  temperature  has  been  modeled  prcviousl)  in  many  v^ays 

Rcieosky,  cl  al.,  used  random  obsuvalums  of  hourly  temperature  daUi  taken  from  o\cr  a  period  of  4  ycar^ 

Rcicosky  independently  tested  and  evaluated  five  different  methods  that  use  daily  maximum  and  minimum 
temperatures  as  input  to  generate  temperatures  for  the  missing  hours  and  concluded  that  the  best  model  ^vas  one 
developed  by  Hoogenboom  and  Huck  (I986).  The  Hwigcnboom/Huck  model  uses  a  cosine  curve  (with  variable 
sunrise)  connecting  maximum  and  minimum  temperatures  to  fit  the  diurnal  temperature  curve.  The  lowest  rocn 
mean  square  error  (RMSE)  achieved  by  this  model  was  1.5  for  one  particular  day  v/hen  the  diurnal  U'cnd  was  smooth 
throughout  the  24-hour  period. 

3.2  The  USAFETAC  Model.  In  USAFETAC’s  slud>,  the  input  temperature  data  dcscril>cd  in  Section  2  was  used 
to  evaluate  three  dilfcrcni  methods  for  filling  a  temperature  curve  given  only  maximum  and  minimum  tcmpcralurcs* 

•  simple  cosine  fit 

•  hypcrlx)iic  tangent  fit 

•  cosine  fit  with  variable  sunrise 


To  lest  each  method,  the  observed  hourly  temperatures  were  compared  to  values  calculated  by  each  method  over  a 
24-hour  interval.  The  results  arc  discussed  in  Section  4. 


3.3  Determination  of  Maximum  and  Minimum  Temperatures.  Maximum  and  minimum  temperatures 
were  dctcmimcd  for  c<ich  station  as  follows:  hourl>  temperature  data  for  the  suition  was  obtained  from 
USAFETAC  s  on-line  relational  database  (DB2).  b  ext.  highest  and  lowest  temperatures  for  each  day  were 
determined  and  .stored  in  a  separate  dataset.  The  dal  .sd  with  the  dail>  maximum  and  minimum  icmpcraluro  wav 
then  used  to  mode!  hourly  temperature  data  as  explained  in  the  following  paragraphs. 


3.4  Simple  Cosine  Fit.  In  ilic  first  method,  a  simple  ctisinc  curve  filling  algorithm  is  applied  to  the  data  As 
shown  in  Figure  I.  the  algorithm  uses  maximum  and  minimum  temperatures  (MAX  I  and  MINI)  for  a  given  day  and 
the  minimum  tcmpcralurc  (MIN2)  for  the  next  day. 


In  the  .simple  cosine  fit  method,  the  minimum  icm|vraiurc  is 
assumed  to  occur  at  O.SOO  local  standard  time  (LST)  and  the 
maximum  at  I4(X)  LST.  For  convctiicncc  in  calculations, 
the  stall  lime  is  set  to  /.cro.  Therefore,  the  ri.se  in 
temperature  from  03(X)  to  1400  l^T  (hours  0  to  9) 
is  described  by  the  portion  of  ilic  cosine 
function  from  A  to  B  in  Figure  1 .  Tlic  fall  from 
I4(X)  10  O.SOO  LST  the  next  day  (hours 
iO  to  2.v)  is  described  by  the  portion 
of  the  cosine  function  from  B  to  (\ 


hours  0  to  9 


24  hours 


Figure  /.  Illustration  of  simple  cosine  fit  to  maximum  and  minimum 
temperatures  (MaXI  and  MINI)  of  one  day  and  minimum 
temperature  (MINI)  of  the  next  day. 
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The  simple  cosine  fii  mclhod  fils  a  o  >lcrm  cosine  Fourier  series  to  each  half  of  the  diunial  curve.  To  bcucr 
understand  this  method  as  well  as  the  one  discussed  in  the  next  section,  it  is  described  in  terms  of  a  linear 
transformation.  An  independent  variable  (jc)  can  be  defined  in  terms  of  the  time  {hr)  .such  that  temperature  become^ 
a  linear  function  of  x  where 


x^cos 


(1) 


for  the  hours  {hr)  from  0  to  9,  and 


X  -  cos 


[ 


{/ir-  10) 


(2) 


for  the  hours  from  10  to  23.  Using  the  equation  for  a  line, 

y  =  fnx+h  (3) 


where 


y  =  dependent  variable  (temperature) 
m  =  slope  of  the  line 

X  =  independent  variable  (a  function  of  time) 
h  =  ordinate  (y  -  intercept), 

and  referring  to  Figure  1 ,  the  slope  for  segment  AB  is  given  by 


-( 


MAX\  -  MIN  \ 

2 


) 


and  for  segment  BC  by 


MAXI -MINI 

2 


It  was  assumed  that  the  >■  -  mlerccpi  value  for  both  segments  was  given  by  half  llic  value  of  the  amplitude  nr 

- “2 -  and  - - ,  respectively.  By  combining  terms,  the  equations  to  takiilaic  lioiirls 

temperatures  for  a  day  become 


V(/ir)=- 


{MAX\-MIN\ 

V  2 


)  •  •  g) 


MAX  1  -I  MIN  1 

+ - T - 


(4( 


where  T  (hr)  is  the  hourly  temperature  for  hours  from  0  to and 

/MAX\-MIN2\  r,.  III  M,U’l  f 

/  (/ir)=  I - '  17  J  * - Y 


MINI 


(5) 


where  /  (hr)  is  the  hourly  temperature  for  hours  from  10  to  2.^.  After  vakulations  were  made,  time  was  tonveried 
back  to  local  standard  time  (LST)  by  adding  5  hours. 
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3.5  Hyperbolic  Tangent  Fit.  in  un  ailcmpl  lo  iinprusc  ^.u^vc  filling,  ihc  andl>.si  iricd  ihc  h>|Krbi)liL  iang».ni 
funtliun.  The  advantage  of  lhi.%  funelion  i.s  dial  il  approaches  the  maMinum  and  minimum  as>nipU)lKall>  m»  dial 
plalcaus,  instead  of  pciiks,  are  reached.  Like  die  cosine,  ihc  maximum  and  minimum  \aucs  arc  both  I.  Using  die 
same  procedures  lo  develop  cqualions  as  in  the  cosine  fii  method,  the  temperature  svas  iransL/nncd  to  a  lin».ar 
function  of  jc,  such  that 

X  =  tank  ^  (b) 

for  die  hours  from  0  to  9,  and 

x=tanh\^ — 73 — ) 

for  the  hours  from  10  lo  23.  Since  die  signs  of  the  slopes  arc  opposite  lho.se  in  ihc  cosine  fil,  the  loUil  cquaiuai  fur 
die  hours  from  0  lo  9  becomes 


MAX\-M!N\  ,  ///^-4.5^  AMA' I -f  A//,V  1 
/  {hr)  = - I  — —  I  + - T - 

and  for  die  hours  from  10  lo  23, 

... ,,  ,  ^fAX  1  -  M!N2  .  /hr  -  !6.5\  MAX  I  +  MINI 

I  (hr)  = - -  j— - •  ({inh  I  - —  j  + - ^ -  (9l 

Again,  lime  is  convened  back  lo  LST  by  adding  5  hours. 

3.6  Cosine  Fit  with  Variable  Sunrise.  In  ihc  temperature  filling  method  that  Rcicosk>.  ci.  aL 
dclcrmmcd  lo  be  ihc  bc.su  sunrise  was  calculalcd  to  the  ncarc-vl  whole  hour  and  that  lime  was  asMgnciI  lo  die 
minimum  lempcralure.  Tlic  maximum  Icmpcralurc  was  assumed  lo  occur  al  !4('0  LST.  To  determine  il  lurih^r 
icmpcralurc  curve  ruling  improvcmcnis  were  possible  in  lhi\  stiid>.  the  cosine  Hi  method  was  UNt.d  with  Minris^. 
t^RlSE)  calculalcd  lo  ihc  nearest  whole  hour.  Sunrise  calculation  is  shown  m  die  Appendix.  The  minimum 
icmpcralurc  vvas  assumed  lo  occur  al  RISE,  ihc  maximum  al  i4(X)  LST.  The  linear  iraiislormaiion  ol  kmperjiur^.  lo 
a  function  of  .t  was 


for  die  hours  from  0  to  (14  -  RISE),  and 

j:  =  ({/.r  -|(  14  -  RISE)  +  11}.  ( l 

lor  die  hours  from  { 14  -  RISC)  +  I  lo  24.  Hour  24  was  calculalcd  to  account  for  when  miiUisc  for  lomorrow  was  an 
hour  later  than  todav.  An>  duplicate  observations  were  removed  when  all  tcnlpcralure^  liad  bccncaleuiaied.  The 
cqualions  lo  calculate  hourly  temperatures  for  a  day  were 

...  MAX  1  -  MIN  I  MAX  I  f  MIN  I 

/  = - 2 - -  4 - -j - 


(12) 


for  the  hours  from  0  to  (14  -  RISE)  where  x  was  obtained  from  Equation  10.  and 


...  MAX  1  -  MINI  MAX  1  +  MINI 
7  = - 2 - -  - - 5 - 

for  the  hours  from  (14  -  RISE)  +  1  to  24,  where  x  was  obtained  from  Equation  11. 
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4.  TEST  RESULTS  AND  PROGRAM. 


4,1  Test  Procedures.  Each  of  tiic  three  temperature  lilling  methods  described  in  Section  3  were  tested  using 
daily  temperatures  for  the  stations  listed  in  Section  2.  Tfie  observed  hourly  temperatures  were  compared  to  the  filled 
temperatures  obtained  from  each  method.  To  meaourc  the  accurac>  of  fit,  the  bias  and  RMSE  were  dctcrmincil. 

The  bias,  or  mean  of  the  error  estimate,  was  determined  from  the  following  equation: 


BIAS  = 


-7/) 

Tv 


(14) 


where  the  difference  between  the  observed  hourly  temperature  and  the  filled  hourly  temperature  (I is  the  error 
of  estimate  and  N  is  the  total  number  of  observations.  A  positive  bias  indicates  that  on  the  average,  the  Idled 
temperatures  were  lower  than  the  observed,  while  a  negative  bias  indicates  filled  temperatures  were  higher  than  the 
observed. 


The  RMSE  was  calculated  from 


RMSE  =  (15) 

The  RMSE  gives  a  measure  of  tlie  si/x  of  the  differences  beiwoin  the  observed  and  filled  temperature  values.  Since 
RMSE  calculation  is  heavily  influenced  by  a  few  large  errors  of  estimate,  RMSE  may  be  larger  than  the  error  fur 
most  of  the  data  points.  To  clarify  this  point,  let’s  suppose  a  hypothetical  day  has  an  error  ol  estimate  of  1 .0  for  each 
hour  except  one,  which  has  an  error  of  estimate  of  8.0.  The  calculated  RMSE  is  about  2.0.  One  large  dilfcrcnce 
between  observed  and  filled  temperature  made  the  other  23  filled  temperatures  appear  to  be  more  in  error  than  they 
actually  were. 

4.2  Results,  The  results  of  hourly  temperature  filling  accuracy  for  all  three  methods  aic  ,hown  in  Figure  ,  units 
arc  dcgrce.s  F. 


Hyperbolic  Cosine  with 

Cosine  Tangent  Sunrise 


station 

BIAS 

RMSE 

BIAS 

RMSE 

BIAS 

RMSE 

Scott  AFB,  IL 

-0.07 

4.23 

-0.07 

4.20 

-0.46 

4.26 

Pease  AFB,  NH 

-0.15 

4.14 

-0.15 

4.10 

-0.30 

4.13 

McChord  AFB,  WA 

-0.15 

3.10 

-0.15 

3. ('6 

-0..50 

3.19 

Luke  AFB,  AZ 

0.0.5 

3.48 

0.06 

3.45 

-0.43 

3.45 

Homestead  AFB,  FL 

-0.23 

2.92 

-0.23 

2.90 

-0.50 

3.04 

Mildenhall  AB,  UK 

0.01 

3.16 

0.01 

3.13 

-0.27 

3. 16 

Ramstcin  AB,  GE 

0.09 

3.36 

0.09 

3. .32 

-0.26 

3.34 

Honolulu,  HI 

-O.fr) 

2.16 

-0.66 

2.14 

-0.92 

2.26 

Figure  2.  Test  results  for  three  temperature  filling  methods.  Units  arc  ®  F. 


As  shown  by  these  results  for  all  monUis  and  all  hours,  the  hyperbolic  tangent  filling  mcilioil  was  the  best  based  on 
RMSE.  It  also  proved  best  when  only  the  typical  summer  pcriotl,  from  June  to  September,  was  u>ed.  This  is  shown 
in  Figure  3. 
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Cosine 

Hyperbolic 

Tangent 

Cosine  with 

Sunrise 

station 

BIAS 

RMSE 

BIAS 

RMSE 

BIAS 

RMSE 

Scott  AFB,  IL 

-0.11 

3.40 

-0.11 

3.40 

-0.53 

3.41 

Pease  AFB,  NH 

-0.26 

3.61 

-0.26 

3.60 

-0.63 

3.63 

McChord  AFB,  WA 

-0.10 

2.70 

-0.10 

2.60 

-0.52 

2.75 

Luke  AFB,  AZ 

0.58 

3.48 

0.58 

3.42 

0.04 

3.45 

Homestead  AFB,  FL 

-0.30 

3.00 

-0.30 

2.97 

-0.53 

3.08 

Mildenhall  AB,  UK 

-0.01 

2.97 

-0.01 

2.94 

-0.36 

2.81 

Ramsiein  AB,  GE 

0.14 

3.10 

0.14 

3.06 

-0.30 

3.07 

Honolulu,  HI 

-0.84 

2.06 

-0.84 

2.05 

-1.07 

2.23 

Figure  3,  Same  as  Figure  2,  except  for  summer  months  from  June  to  September,  Units  arc  ®  F. 


4.3  Median  RMSE  and  Absolute  Extreme  Difference.  Another  way  to  evaluate  aceurae)  of  In  is  to 
examine  the  median  (50lh  percentile)  RMSE  and  absolute  extreme  difference*  Tiicsc  two  statistics  arc  shown  in 
Figure  4  for  each  temperature  filling  method,  for  each  station. 


Cosine 

Hyperbolic 

Tangent 

Cosine  with 

Sunrise 

Station 

Median 

Abs  Ext 

Median 

Abs  Ext 

Median 

Abs  Ext 

Scott  AFB,  IL 

1.81 

39.70 

1.65 

38.76 

1.81 

39.70 

Pease  AFB,  NH 

1.90 

38.70 

1.76 

37.67 

1.90 

38.70 

McChord  AFB,  WA 

1.62 

22.13 

1.51 

22.28 

1.65 

20.41 

Luke  AFB,  AZ 

2.01 

25.15 

1.99 

25.97 

1.90 

26.46 

Homestead  AFB,  FL 

1.59 

23.34 

1.50 

22.57 

1.67 

23.10 

Mildenhall  AB,  UK 

1.68 

22.69 

1.51 

22.06 

1.68 

23.24 

Ramstein  AB,GE 

1.73 

32.40 

1..55 

.’1.67 

1.69 

33.94 

Honolulu,  HI 

1.23 

19.45 

1.17 

19.67 

1.25 

20.00 

Figure  4,  Median  RMSE  and  absolute  extreme  difference  (Abs  Ext)  for  three  temperature  filling 
methods.  Periods  of  record  are  as  in  Figure  2.  Units  are  ®F. 


As  before,  the  hyperbolic  tangent  method  achieved  tlic  lowest  RMSE.  The  median  RMSE  values  eaih  station 
show  that  at  the  5()th  percentile  level,  the  RMSE  ranged  Irom  1.17^  F  in  the  tropics  to  1.99^  F  in  the  mid  latitudes 
At  the  extreme  end,  errors  in  excess  of  an  order  of  magmtiule  greater  than  the  median  eonu-ibulcd  to  raising  ♦he 
overall  RMSE,  as  explained  earlier.  Many  of  the  extreme  errors  arc  caused  by  frontal  passages.  An  ».xamplc  is 
shown  in  Figure  5,  where  the  24-hour  period  containing  the  extreme  for  tlie  hyperbolic  tangent  !ncthod  for  Scott 
AFB  IS  used  to  plot  the  filled  and  observed  hourly  temperatures  for  each  hour.  The  extreme  error  occurred  at  1300 
LSTon  28  January  1977. 
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Temp  (F) 


40 
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0  3  6  9  12  15  18  21  24 

Hour  (LST) 

Observed  Hriy  Temps  Filled  Hriy  Temps 

figure  5.  Plots  of filled  hourly  temperatures  and  observed  hourly  temperatures  for  Scott  AFB,  28  January  1 977, 

In  this  CAumplc,  the  maximum  temperature  for  the  day  actuall>  i)tLurral  at  midnight  and  resulted  in  large  errors  ol 
estimate  (observed  minus  filled  temperature  values)  for  that  da>.  Without  additional  information,  a  computer 
program  cannot  account  for  these  anomalous  cases.  The  effects  of  days  like  the  one  illustrated  are  more  pr .  nuunLcd 
during  seasons  (and  at  locations)  with  frequent  frontal  passages,  but  they  arc  minimiAcd  over  long  periods  of  record. 
An  additional  source  of  variance  will  occur  when  this  method  is  applied  to  the  NOAA  Climatological  Data.  As 
explained  by  McFarland,  Miller,  and  Neale  (1990),  the  reported  minimum  temperature  often  relates  to  the  previous 
day’s  minimum  when  the  observation  time  is  0500  LST. 

4.4  The  Programs.  Based  on  the  favorable  test  results  for  the  hyperbolic  tangent  temperature  filling  method. 
Equations  8  and  9  were  used  as  algorithms  to  write  a  Statistical  Analysis  System  (SAS)  program  for  use  b> 
USAFETAC/ECE  (SAS  is  the  registered  trademark  of  SAS  IndusUics,  Inc,  Cary,  NC).  The  program  (MXMNFILL) 
calculates  values  for  mi.ssing  hourly  temperatures  given  input  of  daily  maximum  and  minimum  temperatures.  To 
account  for  locations  that  close  on  weekends  and  holidays,  another  program  (WEEKEND)  was  developed  to  linearis 
interpolate  daily  maximum  and  minimum  temperatures  when  less  than  .3  consecutive  days  were  missing  For  the  lest 
results  shown  in  Figure  2,  every  6lh  and  7th  day’s  maximum  and  minimum  temperatures  were  set  to  missing.  When 
WEEKEND  was  iisctl  to  replace  the  temperatures,  RMSE  was  raised  by  F  for  the  mid  latitudes  and  subtropics, 
bul  only  slightly  for  the  tropics. 
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5.  SUMMARY. 


5.1  Temperature  Filling  Method.  This  report  describes  the  development  and  testing  of  a  mctliod  for  filling  in 
missing  hourl>  temperatures  at  locations  where  onl>  maximum  and  minimum  dail>  temperatures  are  rcp^)ried.  The 
method  uses  the  hyperbolic  tangent  function  with  the  temperature  iransformul  to  a  linear  (unction  of  the  liour.  The 
advantage  in  using  the  hyperbolic  ta  .ecnl  is  that  temperatures  near  the  maximum  and  minimum  arc  alhjwcd  to 
approach  those  points  gradually. 

5.2  Accuracy.  Tests  of  the  method,  using  hourly  temperature  data  for  all  months  and  all  hours  for  eight 
locations  in  different  climatic  regions,  showed  an  RMSE  from  2^  F  lor  the  tropics  to  4^  F  for  mid-latitudes.  Using 
the  same  data  for  the  typical  summer  period  (June  to  September,  when  frontal  passage  is  Ic.ss  frequent),  RMSE 
values  were  about  10  percent  lower.  If  this  method  is  applied  to  stations  that  only  report  maximum  and  minimum 
temperatures,  accuracy  will  decrease  further  due  lo  variances  associated  with  this  t>pc  of  reporting  station. 

5.3  The  Programs.  USAFETAC/DNO  used  the  hourly  temperature  algorithms  developed  lor  the  h>pcrlx)lie 
tangent  temperature  filling  method  to  write  two  SAS  programs  for  USAFETAC/ECE.  One  program  (MXMNFILL) 
generates  a  dataset  of  estimated  hourly  temperatures  given  input  of  daily  maximum  and  minimum  tcmperaiup.". 
The  other  (WEEKEND)  accounts  for  locations  that  may  not  be  open  7  days  a  week  b>  providing  micrpolaied  dailv 
maximum  and  minimum  temperatures  for  cases  in  which  the  data  gaps  are  less  than  ^  days. 
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APPENDIX 


Sunrise  Calculations 


The  sun’s  declination  angle,  DECLIN,  is  obtained  from  Cox  (1987): 

DECUN  =  -23.45  cos  \{JULSEC  +  777,600)  (fRTORAD)] 

where  JULSEC  is  the  number  of  seconds  from  I  January  I960  and  777, 6(X)  is  the  number  of  seconds  from  the  inter 
solstice  to  1  Janiuiry;  YRl  GRAD  =  ^3  ^24  >1  60()  account  the  fact  that  die  earth  revolves 

around  the  sun  (2K)  in  1  year.  The  solar  zenith  angle  (Z)  can  be  calculated  from: 

cos  (Z)  =  sin  {LA'n  stn{DECUN)  +  cosiLAT)  cosiDECUN)  cosiHRANG ) 

where  LAT  is  the  latitude  of  the  site,  and  HRANG  is  the  hour  angle  of  the  sun.  At  sunrise,  the  zenith  angle  is  -90® 
or  cos(Z)  =  0.  Solving  the  above  equation  for  HRANG  yields 

COSilRANG  =  -  (TANDEC  •  TANLAl ) 

where  the  tangent  of  the  declination  angle, 

TANDEC  =  TAN  (dECLIN  •  ■^) 

and  the  uuigeni  of  the  latitude  in  decimal  form  (LATDEC)  is  given  by 

TANLAT  =  TAN  (^LATDEC  •  -j^) 

When  COS  HRANG  is  grcalcr  than  -1  and  less  than  or  equal  to  1 ,  the  sun  would  never  rise.  For  our  purposes 
however,  we  a.ssunic  sunrise  is  at  noon  so  that  a  sunri.se  to  the  14(X)  LST  portion  of  the  curve  can  be  calculated 

HRANG  =ARCOS  (COSilRANG) 


Otherwise,  HRANG  =  0.  The  sunrise  in  Uxal  time  is  then  obtained  from 


RISE=  12- 


HRANG 


where  the  sun  is  vertical  at  1200  LST  and  1  hour  of  time  is  equivalent  U)  15"  longitude. 


II 
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GLOSSARY 


Abs  Exl  absolute  extreme  difference 

CDH  C(K)ling  degree  hours 

DB2  US  AFETAC’s  on-line  relational  database 

DNO  the  Operations  Applications  Development  Section  of  USAFETAC’s  Aerospace  Sciences  Branc  h 

ECE  the  Engineering  Meteorology  Section  of  USAFETAC’s  Environmental  Applications  Branch 

F  Fahrenheit 

hr  hour 

LST  local  standard  time 

MAX  I  maximum  temperature  of  1  day 

MINI  minimum  temperature  of  1  day 

MIN2  minimum  temperalurc  of  the  next  day 

FOR  period  of  record 

RISE  sunrise 

RMSE  root  mean  square  error 
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